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Abstract 

The unsteady, incompressible Navier-Stokes equa- 
tions are used for the direct numerical simulation 
(DNS) of spatially evolving disturbances in a three- 
dimensional (3-D) attachment-line boundary layer. 
Two-dimensional (2-D) disturbances are introduced ei- 
ther by forcing at the inflow or by harmonic-source gen- 
erators at the wall; 3-D disturbances are introduced 
by harmonic-source generators at the wall. The DNS 
results are in good agreement with both 2-D non- 
parallel theory (for small-amplitude disturbances) and 
weakly nonlinear theory (for finite-amplitude distur- 
bances), which validates the two theories. The 2-D DNS 
results indicate that nonlinear disturbance growth oc- 
curs near branch II of the neutral stability curve; how- 
ever, steady suction can be used to stabilize this dis- 
turbance growth. For 3-D instabilities that are gen- 
erated off the attachment line, spreading both toward 
and away from the attachment line causes energy trans- 
fer to the attachment-line and downstream instabilities; 
suction stabilizes these instabilities. Furthermore, 3-D 
instabilities are more stable than 2-D or quasi-2-D in- 
stabilities. 

1. Introduction 

Many instability mechanisms can occur that cause 
the breakdown from laminar to turbulent flow on swept 
wings; however, this discussion will focus on only those 
disturbances that evolve near the attachment-line re- 
gion (near the leading edge). Turbulent contamination, 
which results from turbulence at a fuselage- wing junc- 
ture, can travel out over the wing and cause laminar 
flow on the wing to become turbulent. To prevent this 
contamination, devices such as the Caster bump 1 or suc- 
tion (Pfenninger 2 ), implemented near the wing root, can 
halt the turbulent attachment-line boundary layer from 
sweeping out over the entire wing. 

Although the problem of turbulent flow that orig- 
inates from the fuselage-wing juncture can be avoided 
with a device such as the Caster bump, a Reynolds num- 
ber must exist beyond which disturbances generated by 
surface imperfections or particulates on the wing, cou- 
pled with noise, eventually cause transition. If we 
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assume that the initiated disturbances are sufficiently 
small, then hydrodynamic stability theory could poten- 
tially be used to predict the spatial amplification and 
decay of the disturbances along the attachment line. 
Caster 3 first examined this small-amplitude disturbance 
problem by using acoustic excitation along the attach- 
ment line of a swept-cylinder model. Caster generated 
sine waves with various frequencies that were detected 
in the flow by a hot-film gauge on the attachment line. 
He noted that the recorded oscillations had preferred 
frequency bands that changed with tunnel speed and 
that this behavior was reminiscent of traveling-wave in- 
stabilities. From his measurements, he concluded that 
the small-amplitude disturbances in an attachment-line 
boundary layer were stable for momentum-thickness 
Reynolds numbers Rg below 170 (the critical Reynolds 
number was outside the experimental range). This 
value for the critical Reynolds number is close to the 
theoretical value of 200, which is obtained by assum- 
ing a two-dimensional (2-D) attachment-line bound- 
ary layer. Later, Cumpsty and Head 4 experimentally 
studied large-amplitude disturbances and turbulent flow 
along the attachment line of a swept-wing model. With- 
out tripping the boundary-layer instabilities, they ob- 
served that laminar flow is stable to small-amplitude 
disturbances up to Rg ~ 245 (which corresponds to the 
top speed of the tunnel). They note that this observa- 
tion does not contradict the theoretical value; an accu- 
rate theoretical value would need to account for three- 
dimensional (3-D) effects. At the same time, Pfenninger 
and Bacon 5 used a wing swept to 45° to experimentally 
study the attachment-line instabilities in a wind tunnel 
that could reach speeds large enough to obtain unstable 
disturbances. With hot wires, they observed regular si- 
nusoidal oscillations with frequencies comparable to the 
most unstable 2-D modes of theory; these modes caused 
transition at about Rg ~ 240. A continued interest in 
the transition initiated near the attachment line of swept 
wings led Poll 6,7 to perform additional experiments with 
the swept circular model of Cumpsty and Head 4 . Like 
Pfenninger and Bacon 5 , Poll observed disturbances that 
amplified along the attachment line. He noted that no 
unstable modes were observed below Rg = 230. 

With nonparallel stability theory, Hall et al. 8 stud- 
ied the linear stability of the attachment-line boundary- 
layer flow called swept Hiemenz flow, which is sketched 
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in Fig. 1. This 3-D base flow is an exact solution of 
the Navier-Stokes equations; hence, its use is advanta- 
geous in stability analyses. By assuming 2-D instability 
modes, Hall et al. 8 demonstrated that the attachment- 
line boundary layer can theoretically be stabilized with 
small amounts of suction. Hall and Malik 9 extended this 
theory to account for weakly nonlinear stability. They 
note that subcritical instability is observed at wave num- 
bers that correspond to the upper branch of the neu- 
tral curve. The neutral curve and the subcritical insta- 
bility region of Hall and Malik 9 are shown in Fig. 2 
with the experimental results of Pfenninger and Bacon 5 
and Poll. 5 -' Consistent with the Pfenninger and Bacon 5 
experimental results, large-amplitude disturbances be- 
came unstable before the linear critical point. 
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Fig. 1. Sketch of attachment-line region of swept 
Hiemenz flow. 

Hall and Seddougui 10 studied oblique waves and 
their interaction in attachment-line flow at the large 
Reynolds number limit. They nop that close to the 
attachment line a small band of destabilized oblique 
modes appears, interacts with the 2-D mode, and causes 
a breakdown of the 2-D mode. Further, they nofe that 
oblique modes become less important further away from 
the attachment line and that low-frequency modes be- 
come the dominant mechanism. (Perhaps stationary 
crossflow modes arise near and dominate away from the 
attachment line.) 

Finally, Spalart 11 used a direct numerical simula- 
tion (DNS) approach, based on the fringe method, to 
study the leading-edge contamination problem. Small- 
amplitude disturbances were initialized with white 
noise. A Reynolds number range with a width of 50 
was selected. At the lower Reynolds number, all distur- 
bances decayed; at the higher number, at least one mode 
was amplified. The neutral curve predicted by Hall 
et al. 8 fell within the Reynolds number range used by 


Spalart. These results show qualitative agreement with 
the linear theory of Hall et al. 8 Furthermore, Spalart, 11 
demonstrated that simple Hiemenz flow (i.e. , R = 0) is 
both linearly and nonlinearly stable. 
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Fig. 2. Neutral curve, experimental regions of instabil- 
ity growth, and theoretical region of subcritical growth 
in attachment-line boundary layer. 

For the present study, a 2-D spatial DNS code de- 
scribed by Joslin et al. 12,13 and a newly developed 3- 
D spatial DNS code described in this paper are used 
to study the linear and nonlinear instabilities that ini- 
tiate and develop along the attachment line of swept 
Hiemenz flow. The region near branch II of the neutral 
curve is investigated with DNS to verify the nonparal- 
lel theory of Hall et al. 8 for infinitesimal disturbances. 
Subcritical instability growth in the region shown in 
Fig. 2 was predicted with the weakly nonlinear the- 
ory of Hall and Malik; 9 however, recent DNS results 
by both Jimenez et al. 14 and Theofilis 15 failed to find 
this subcritical growth. Jimenez et al. 14 contend that 
this subcritical growth region does not- exist. Here, an 
independent study is conducted to determine if nonlin- 
early growing disturbances can be found. Furthermore^ 
because the true physical problem involves 3-D distur- 
bances that may be on or off the attachment line, 3-D 
simulations are conducted to determine the importance 
of 3-D disturbances near the attachment line. Steady 
suction is used to control instability growth. 
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2. Problem Formulation 

For the problem at hand, the velocities u = (u, v, w) 
and the pressure p are solutions of the incompressible, 
unsteady Navier-Stokes equations. The instantaneous 
velocities u and the pressure p may be decomposed into 
base and disturbance components as 

{£, P}(x, t) = { U_, P}{x) + {u, p}(x, t) (1) 

The base flow is given by the velocities U_ = ( U , V, W ) 
and the pressure P ; the disturbance component is given 
by the velocities u = ( u,v,w ) and the pressure p. A 
Cartesian coordinate system x_ = ( x,y,z ) is used, in 
which x is aligned with the attachment line, y is wall 
normal, and z corresponds to the direction of flow ac- 
celeration away from the attachment line. 

Originally described by Hall et al. 8 , the base flow 
referred to as a swept Hiemenz flow is an exact solution 
to the incompressible 3-D Navier-Stokes equations. The 
velocities at the wall are 

U = W = 0 and V = V 0 at y = 0 (2) 

Away from the wall, the velocities are 

U — >■ U o and W —>■ W 0 — as y — >■ oo (3) 

±j 

where U 0 , V 0 , W 0 are independent velocity scales and L 
is the length scale in the flow-acceleration direction z. 
Shown in Fig. 1, the fluid comes straight down toward 
the wall; it turns away from the attachment line into 
the ±z directions to form a boundary layer. In the x 
direction, the flow is uniform. In the absence of sweep, 
U 0 is equal to 0, and the flow reduces to the 2-D stagna- 
tion flow first described by Hiemenz. 16 A boundary-layer 
thickness is defined in the t/z-plane as 6 = \J vL/W 0 \ a 
Reynolds number, as R = IJ 0 b /v = 2.475/?#; and a tran- 
spiration constant, as k = V 0 \/LjvW 0 , where k = 0 for 
the zero-suction case. If the attachment line is assumed 
to be infinitely long, the velocities become functions of 
z and y only, and an exact solution of the Navier-Stokes 
equations can be found. 

In accordance with Hall et ah, 8 the equations for 
the base flow are 

W+V Y = 0 (4) 

Vyyy + by — VVyy — 1 = 0 (5) 

Uyy - VU Y = 0 ( 6 ) 

where {X, Y, Z} = {x,y,z}/6, the subscript Y denotes 
derivatives with respect to Y , and the hats refer to sim- 
ilarity variables. The boundary conditions are given by 

V Y = 0, V = k, and U = 0 at Y = 0 (7) 


Vy — * — 1 and U —> 1 as Y —> oo (8) 

If the solutions of equations (4)-(8) are nondimensional- 
ized with respect to the attachment-line velocity U 0 , the 
boundary-layer thickness 6, and the kinematic viscosity 
v, then the base flow is 

U(Y) = U(Y), V(Y) = 1 1 >(Y), W(Y,Z) = ^W(Y) 

(9) 

For the disturbance portion of equation (1), the 3- 
D incompressible Navier-Stokes equations are solved in 
disturbance form as 

du 1 

^ + (u-V)u + (U-V)u + (u-V)U = -A p+-X 2 u (10) 
ut K 

with the continuity equation 

V-« = 0 (11) 

and boundary conditions 

u = 0 at Y = 0 and u — >■ 0 as Y — >■ oo (12) 

Disturbances are forced either at the inflow or 
by harmonic-source generators, which use suction and 
blowing at the wall and are assumed to decay to zero in 
the far Held. At the inflow, solutions of the base flow are 
forced, and the buffer-domain technique is employed as 
the outflow condition. 

3. Numerical Methods of Solution 

In the attachment-line (A) direction, discretization 
is accomplished with fourth-order central finite differ- 
ences for the pressure equation and sixth-order compact 
differences for the momentum equations in the interior of 
the computational domain. At the boundary and near- 
boundary nodes, fourth-order forward and backward dif- 
ferences are used. The discretization yields a pentadiag- 
onal system for the finite-difference scheme and a tridi- 
agonal system for the compact-difference scheme. The 
approximations can be solved efficiently by appropriate 
backward and forward substitutions. 

In both the wall-normal (Y) and flow-acceleration 
(Z) directions, Chebyshev series are used to approx- 
imate the disturbances at Gauss-Lobatto collocation 
points. A Chebyshev series is used in the wall-normal 
direction to provide good resolution in the high-gradient 
regions near the boundaries. Furthermore, the use of as 
few grid points as possible results in significant com- 
putational cost savings. In particular, the use of the 
Chebyshev series enables an efficient pressure solver. 
Because this series and its associated spectral operators 
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are defined on [-1, 1] and because the physical prob- 
lem of interest has a truncated domain [0,j/max] and 
[— Zmax, Zmax], transformations are employed. Further- 
more, stretching functions are used to cluster the grid 
near both the wall and the attachment line. For fur- 
ther details on the properties and the use of spectral 
methods, refer to Canuto et al. 17 

For time marching, a time-splitting procedure was 
used with implicit Crank-Nicolson differencing for nor- 
mal diffusion terms; an explicit three-stage Runge-Kutta 
(RK) method was used for the remaining terms. The 
pressure is omitted from the momentum equations (10) 
for the fractional RK stage, which leads to 

rhi * 1 

-= +(m*-V)m* + (17-V)m*+(m*-V) 17 = -VV (13a) 
ut K 

h m h m 

l£=u BC + h? (l + -^r T )Vp™--f rT Vp™- 1 
ht "t 

(m 

where h™ are time-step sizes in the RK scheme, time is 
advanced from u m to the intermediate disturbance ve- 
locities u * , and u* are intermediate boundary conditions 
described by Joslin et al. 12 to reduce the time-marching 
slip velocities that may arise. For the boundary con- 
ditions, u BC = 0 for a rigid wall and u BC = u 0 for an 
inflow condition or a wall-slot condition evaluated at the 
appropriate time in the RK stage. 

A full RK stage is completed by advancing the so- 
lution in time from u* to u m+1 by 

Fhi m+1 

1= = -W" +1 (14a) 

V-« m+1 =0 (146) 

If the divergence of equation (14a) is taken and zero 
divergence of the flow Held is imposed at each RK stage, 
then a pressure equation is obtained: 

V 2 p m+1 = ^(V-«*) (15) 

which is subject to homogeneous Neumann boundary 
conditions. This boundary condition is justified in 
the context of a time-splitting scheme as discussed by 
Streett and Hussaini. 18 

The solution is determined on a staggered grid. 
The intermediate RK velocities («*) are determined by 
solving equation (13) on Gauss-Lobatto points. The 
pressure (p m+1 ) is found by solving equation (15) on 
Gauss points and is then spectrally interpolated onto 
Gauss-Lobatto points. Then, the full RK stage veloci- 
ties (w m+1 ) are obtained from equation (14) on Gauss- 
Lobatto points. The above system is solved three con- 
secutive times to obtain full time-step velocities. The 


three-stage RK time steps given by Williamson 19 are 
{h], h 2 t , h 3 } = {1/3, 5/12, 1/4 }/i t , where the sum of the 
three RK time stages equals the full time step (/i t ). 

To satisfy global mass conservation, an influence- 
matrix method is employed and is described in some de- 
tail by Streett and Hussaini, 18 Danabasoglu, Biringen, 
and Streett; 20,21 and Joslin et al. 12,13 For boundary- 
layer flow, four Poisson-Dirichlet problems are solved 
for the discrete mode that corresponds to the zero eigen- 
value of the system; single Poisson-Neumann problems 
are solved for all other modes. To efficiently solve the 
resulting Poisson problem given by equation (15), the 
tensor-product method of Lynch et al. 22 is used. The 
discretized form of equation (15) is 

(l x (giI(giI + I(giL y (giI + I(giI(gi L z ^jp = R (16) 

where p is the desired pressure solution; R is the right 
side of equation (15); I is the identity matrix; L x is the 
attachment-line-directed central finite-difference opera- 
tor; L y and L z are the wall-normal-directed and flow- 
acceleration-directed spectral operators; and (gi infers a 
tensor product. By decomposing the operators L y and 
L z into their respective eigenvalues and eigenvectors, we 
find 

L y = QAyQ- 1 and L z = SA Z S~ 1 (17) 

where Q and S are the eigenvectors of L y and L z , Q -1 
and S ' -1 are inverse matrices of Q and S, and Ay and 
A z are the eigenvalues of L y and L z . The solution pro- 
cedure reduces to the following sequence of operations 
to determine the pressure p: 

p* = (I <X> Q -1 <X> S~ 1 )R 

pt = (L x (g I (g I + I (g A y (g I + I (g I (g A z )~ 1 p* 
p = (I <X> Q <X> S)p^ (18) 

Because the number of grid points in the attachment- 
line direction is typically an order of magnitude larger 
than the wall-normal and flow-acceleration directions, 
the operator L x is much larger than both L y and L z . Be- 
cause L x is large and has a sparse pentadiagonal struc- 
ture and because A y and A z influence the diagonal only, 
an LU decomposition is performed for the second stage 
of equation (18) once, and forward and backward solves 
are performed for each time step of the simulation. The 
first and third steps of the pressure solve for equation 
(18) involve matrix multiplications. 

To obtain the attachment-line-directed operator 
L x , central finite differences are used. To find the wall- 
normal L y and flow-acceleration L z operators, the fol- 
lowing matrix operations are required: 

Ly = I% L DybyI% L UUd L Z = I^D Z D Z 1^ {IS) 
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where D y is a spectral, wall-normal derivative opera- 
tor for the stretched grid; D z is the spectral, derivative 
operator that is grid clustered in the attachment-line re- 
gion; and D y and D z are the derivative operators with 
the first and last rows set to 0. The interpolation matrix 
I° L operates on variables at Gauss-Lobatto points and 
transforms them to Gauss points; the interpolation ma- 
trix I° L performs the inverse operation. The spectral 
operators are described in detail by Canuto et al. 17 and 
Joslin et al. 4 * * * * * * * * 13 

The operators {L x , L y , L z } , the eigenvalue matri- 
ces {Ay, A^}, the eigenvector matrices {Q, Q -1 , S, S' -1 }, 
and the influence matrix are all mesh-dependent matri- 
ces and must be calculated only once. 

The buffer-domain technique introduced by Streett 
and Macaraeg 23 is used for the outflow condition. As 
shown by Joslin et al. 12 for the flat-plate boundary- 
layer problem, a buffer length of three disturbance wave- 
lengths is adequate for traveling waves. The distur- 
bances are assumed to be from the discrete spectrum, 
which exponentially decays with distance from the wall. 
Both at the wall and in the far Held, homogeneous 
Dirichlet conditions are imposed. For 3-D simulations, 
both homogeneous Dirichlet and Neumann conditions 
have been used in the flow-accelerated direction. With 
either condition, the disturbance will develop in the 
same manner along the attachment line, provided that 
the boundaries are sufficiently far from the attachment- 
line region. The base flow is used for the inflow bound- 
ary condition. Finally, disturbances are forced by un- 
steady suction and blowing of the wall-normal veloc- 
ity component through the wall (harmonic-source gen- 
erators). Although the disturbances may be generated 
by random frequency input, the disturbances of interest 
here are forced with known frequencies. 

4. Two-Dimensional Simulation Results 

The simulations are performed on a grid of 661 

points (~ 60 points per wavelength) along the attach- 

ment line and 81 points in the wall-normal direction. 

The far-field boundary is located at 50 6 from the wall, 

and the computational length along the attachment line 

is 216.56 6. For the time-marching scheme, the distur- 
bance wavelength was divided into 320 time steps per 
period for small-amplitude disturbances and into 2560 
time steps for large-amplitude disturbances (stability 

considerations). The total Cray Y-MP time for a sim- 
ulation with a single processor was 1.5 hrs for small- 

amplitude disturbances and 11 hrs for large-amplitude 

disturbances. 

4.1 Subcritical Growth 

Recall that Hall and Malik 9 developed a weakly 


nonlinear theory to predict the growth and decay of in- 
stabilities in the attachment-line boundary-layer flow. 
With swept Hiemenz flow, they discovered regions where 
2-D disturbances encountered subcritical growth; how- 
ever, the linear theory predicted stable modes. Later, 
both Jimenez et al. 14 and Theofilis 15 discounted this dis- 
covery with recently developed 2-D DNS codes. To re- 
solve this discrepancy, the results from the present study 
of attachment-line boundary-layer disturbances in swept 
Hiemenz flow (computed with a well-tested spatial DNS 
code) are compared with the previous studies of Hall et 
al. 8 , Hall and Malik 9 , Jimenez et al. 14 , and Theofilis. 15 
In addition, the effects of suction on unstable modes are 
documented. 

Disturbances for the first simulations are forced at 
the computational inflow with an amplitude of A = 
0.001 percent (i.e. , some small amplitude), to compare 
the DNS results with linear stability theory (LST). The 
Reynolds number R = 570 and the frequency u> = 0.1249 
correspond to the region of subcritical growth found 
by Hall and Malik 9 , in which disturbances are linearly 
stable. Disturbances that evolve in both a base flow 
that complements the quasi-parallel LST assumptions 
( V = 0) and the full, swept Hiemenz flow are computed 
with DNS. The computed disturbance decay rate and 
the wavelength in the quasi-parallel flow agree with LST. 
The disturbance that propagates in the complete swept 
Hiemenz flow closely retains the wavelength predicted 
by LST, but decays at a slower rate than that predicted 
by LST. This change in decay rate is consistent with 
the theory of Hall et al. 8 From this comparison, we find 
that the wall-normal velocity (V) terms in the stability 
equations have a destabilizing effect on the disturbance. 

The fundamental wave, the mean-flow distortion, 
and the harmonics from a simulation forced at the inflow 
with a large amplitude of A = 12 percent, a Reynolds 
number of R = 570, and a frequency of u> = 0.1249 
are shown in Fig. 3. After a transient region of ad- 
justment, the fundamental wave encounters subcritical 
growth, which is in agreement with the weakly nonlin- 
ear theory of Hall and Malik. 9 Because DNS involves 
many parameters (e.g., grid, time-step size, order of 
methods), the potential for errors was high in the work 
of Jimenez et al.; 14 erroneous results may have misled 
them to refute the concept of subcritical growth. For 
this flow, instantaneous and mean attachment-line and 
wall-normal velocity profiles at various attachment-line 
locations are shown in Fig. 4. The profiles indicate 
that time-dependent distortions to the base flow are ob- 
served, but the mean flow (U + u 0 ), which consists of 
the base flow and the mean-flow distortion components, 
shows little deviation from the base-flow solution. How- 
ever, both the time-dependent and mean wall-normal 
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profiles undergo distortions because of the disturbance. 
Finally, Fig. 5 shows the wall-normal component of the 
base flow that corresponds to R = 570 and R = 670. 
A comparison of these base-flow profiles with the mean 
flow of Fig. 4 shows that a large-anrplitude disturbance 
produces a distortion to the base flow, which causes an 
effective increase in the base Reynolds number. The 
increase in Reynolds number alone does not account 
for the growing mode (based on linear stability analysis 
with the same frequency). However, we surmise that 
(similar to nonparallel effects) nonlinear disturbances 
broaden the neutral curve toward higher frequencies in 
the critical region. 



Fig. 3. Nonlinear subcritical disturbance growth in 
attachment-line boundary layer at R = 570 and u> = 


0.1249. 
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Fig. 4. Instantaneous and mean velocity profiles of non- 
linear, subcritically growing disturbance in attachment- 
line boundary layer at R = 570 and u> = 0.1249. 
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Fig. 5. Wall-normal component of base flows that cor- 
responds to R = 570 and R = 670. 

To control the subcritical growth of disturbances, 
various levels of suction are employed. Although Hall 
and Malik 9 note that suction makes the flow more sus- 
ceptible to subcritical instability growth, Fig. 6 shows 
that this subcritical disturbance growth can be con- 
trolled by using small levels of suction. If the 2-D DNS 
results mimic the actual 3-D behavior of the flow, then 
large-anrplitude disturbances generated on the attach- 
ment line can be controlled with suction. 



Fig. 6. Control of nonlinear subcritical disturbance 
growth in attachment-line boundary layer at R = 570 
and u> = 0.1249 with suction. 
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4.2 Linear fc Nonlinear Instability at Upper Branch 

Stability in the region near branch II of the neutral 
curve is studied. The disturbances are forced by suction 
and blowing at the wall with an amplitude of A = 0.001 
percent: The growth and the decay of various frequency 
waves at R = 684.2 are shown in Fig. 7. Approximate 
growth rates are listed in the figure with a superimposed 
neutral solution (horizontal line). These results are in 
agreement with the neutral curve that was predicted 
theoretically by Hall et. a.l. 8 The results also agree with 
a recent linear DNS study by Theofilis. 24 
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Fig. 7. Instability growth and decay near neutral curve 
of attachment-line boundary layer at R = 684.2. 


The effect of both steady suction and steady blow- 
ing on linear instability growth in this region is docu- 
mented. The amplification of a growing wave at R = 
684.2 and u> = 0.1150 is shown with the growth or de- 
cay of similar waves that evolve in swept Hiemenz flow 
in the presence of steady suction or blowing. The results 
indicate that suction stabilizes and blowing significantly 
destabilizes the disturbance. The effects of suction and 
blowing on disturbances computed by DNS are in agree- 
ment. with the theory of Hall et al. 8 for small-a.mplitude 
disturbances. 

To determine if nonlinear disturbance growth or 
equilibrium states can be found above branch II of the 
neutral stability curve and to ensure that, the subcrit.i- 
ca.l growth (obtained both by Hall and Malik 9 and the 
DNS shown in Fig. 3) did not. artificially result, from 
the disturbance forcing at. the inflow boundary, a. se- 
quence of simulations was performed with harmonic- 
source generators at. the wall. The disturbances are 
forced by suction and blowing at. a. decaying mode that, 
corresponds to R = 684.2 and u> = 0.1249. The ini- 


tial amplitudes of the disturbances for each simulation 
were incrementally increased until the otherwise linearly 
decaying mode became amplified because of the nonlin- 
ear forcing. With steady suction, the results further 
demonstrate that, small amounts of suction can be used 
to stabilize these disturbances. To study the subcrit.- 
ica.l growth phenomenon, Theofilis 15 used suction and 
blowing to initiate disturbances in the simulations. It. is 
clear from the present, results that. Theofilis 15 failed to 
obtain subcrit.ica.l growth because the suction and blow- 
ing amplitudes were too small to generate a. subcrit.ica.lly 
growing disturbance. 

The simulations are repeated for the mode that, cor- 
responds to R = 684.2 and u> = 0.1230, which is closer 
to branch II of the neutral curve. The results indi- 
cate that, nonlinear growth can be obtained from dis- 
turbances with smaller initial amplitudes. As expected, 
smaller amounts of suction are required to stabilize the 
disturbances. This series of simulations was conducted 
to demonstrate nonlinear growth close to branch II of 
the neutral curve; these results create a. foundation for 
3-D nonlinear simulations. Because significant, compu- 
tational cost, is involved with 3-D simulations and be- 
cause the forcing of disturbances with very large ini- 
tial amplitudes requires more grid resolution, these 2-D 
simulations have established that, forcing smaller ampli- 
tude disturbances close to branch II will yield nonlin- 
ea.rly growing modes. Furthermore, smaller amounts of 
suction can be used to stabilize these growing modes. 

4.3 Linear fc Nonlinear Instability at. Lower Branch 

Finally, the region near branch I of the neutral 
curve is studied. The evolution of disturbances gen- 
erated by suction and blowing at. R = 684.2 and u> = 
0.0840 is shown in Fig. 8. For sma.ll-amplit.ude distur- 
bances, the theory of Hall et. al. 8 give slightly different, 
results compared with both the present. DNS results and 
the DNS results of Theofilis. 24 Because the relative er- 
rors between the theory and DNS frequencies are only 
1 to 2 percent., the predit.ions of Hall et. al. 8 for the in- 
stability characteristics are correct, for practical applica- 
tions. Various large-amplitude disturbances were forced 
by suction and blowing; no subcrit.ica.l growth was de- 
tected near branch I of the neutral curve. However, 
supercritical la.rge-a.mplit.ude disturbances approach a. 
state that, agrees with Hall and Malik. 9 Furthermore 
this equilibriumlike state agrees with the experimental 
observations of Pfenninger and Bacon. 5 

5. Three-Dimensional Simulation Results 

The 3-D simulations are performed on a. grid of 661 
points along the attachment, line, 81 points in the wall- 
normal direction, and 25 points in the flow-a.ccelera.t.ion 
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direction. The far-held boundary is located at 50(5 from 
the wall, the computational length along the attachment 
line is 216. 56(5, and the how-acceleration boundaries art 
located ±1005 from the attachment line. For the time- 
marching scheme, the disturbance wavelength was di- 
vided into 320 time steps per period. The total Cray 
Y-MP time for a simulation with a single processor was 
approximately 25 hrs (with a single processor). 
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Fig. 8. Linear decay and nonlinear equilibrium in 
attachment-line boundary layer at R = 684.2 and 
l o = 0.0840. (Disturbances normalized by initial am- 
plitudes.) 

5.1 Region of Disturbance Decay 

The nonparallel theory of Hall et al. 8 outlined 
the stable and unstable regions for infinitesimal distur- 
bances. In a segment of the subcritical region, large- 
amplitude disturbances were found by Hall and Malik 9 
to exhibit nonlinear amplification. The 2-D DNS re- 
sults confirmed this subcritical growth phenomenon. In 
this section, the Reynolds number R = 570 and the 
frequency u> = 0.1249, which are parameters in the sub- 
critical region, are used in the study of the evolution 
of small-amplitude 3-D disturbances. The results are 
compared with LST and the 2-D DNS results. 

To compare with the 2-D results, a quasi-2-D dis- 
turbance is initiated in the 3-D flow. At best, this dis- 
turbance is an approximation to a true 2-D instability. 
To generate this 2-D disturbance, a harmonic source is 
used that is elongated (—44.2 < Z < 44.2) in the flow- 
acceleration direction. This disturbance-forcing method 
is comparable to using a vibrating ribbon to generate 2- 
D disturbances for wind-tunnel experiments. The qual- 
itative features of a disturbance generated by the har- 
monic source with a small amplitude (e.g., A = 0.001 



percent) are shown in Fig. 9. The disturbance evolu- 
tion is viewed from above and along the attachment line. 
The wave travels along the attachment line without sig- 
nificant 3-D features. However, because the base flow is 
accelerating away from the attachment line (in the ±Z 
directions), wave spreading occurs with distance from 
the harmonic source, and the raff of spreading increases 
with distance along the attachment line. 


x 1 CD — 
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Fig. 9. Side and top view of 3-D traveling wave in 
attachment-line boundary layer. 


The quasi-2-D simulation results for both a quasi- 
parallel base flow (i.e. , V = 0) and the full swept 
Hiemenz flow were compared with linear stability the- 
ory. The amplitude, decay rate, and wavelength of dis- 
turbances simulated with the quasi-parallel flow are in 
excellent quantitative agreement with the 2-D LST re- 
sults. This agreement suggests that in this parameter 
region the elongated harmonic source can approximate 
a 2-D disturbance on the attachment line. Similar to 
the 2-D DNS comparison, the full swept Hiemenz base 
flow destabilizes disturbances due to the inclusion of the 
V velocity component. This destabilizing feature is con- 
sistent with theory. 

To further demonstrate the 2-D nature of the dis- 
turbance generated with the elongated harmonic source, 
Fig. 10 shows the attachment-line results compared 
with results at a distances 135 and 355 off the attach- 
ment line. The evolution patterns are identical out to 
near 355, where small deviations are observed. This im- 
plies that the elongated harmonic source is generating 
primarily 2-D waves and that the attachment-line ve- 
locity component is dominant (i.e., the amplitude of the 
w velocity component of the disturbance is too small 
to modify the dominant u component). The u and w 
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velocity profiles at Z = 135 and 355 reveal only small 
differences with u velocity components and large differ- 
ences with w velocity components. The w velocity is an 
order of magnitude smaller than the u velocity, which is 
the reason for the good agreement between the u veloc- 
ity on the attachment line with the same components 
off the attachment line. Furthermore, although no sym- 
metry assumption is made, flow symmetry about the at- 
tachment line is realized with this particular harmonic- 
source generator. 
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Fig. 10. Flow-acceleration variation of simulated 2-D 
disturbance evolution in 3-D attachment-line base flow 
for R = 570 and w = 0.1249. 
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the 2-D simulations, amplitude information is overpre- 
dicted, and growth-rate information is underpredicted. 
These results suggest that much larger disturbances will 
be required to generate subcritical disturbance growth 
in the 3-D flow (if subcritical growth is possible in the 
3-D flow). 



X 



x 

Fig. 11. Comparison of 2-D and 3-D disturbance evolu- 
tions in 3-D attachment-line boundary layer for R = 570 
and lo = 0.1249. 


In Figs. 11 and 12, 3-D simulation results on the 
attachment line are compared with previous 2-D DNS 
results. Figure 11 clearly shows a significant amplitude 
disparity between the 2-D and 3-D results. Because 
the 3-D simulations contain a flow-acceleration veloc- 
ity component w, an additional degree of freedom is 
available to disperse (or absorb) energy. Hence, the 
harmonic-source generator forces less energy into the 
attachment-line velocity component u. The 2-D and 3-D 
results (normalized by the 2-D maximum of the u veloc- 
ity) are also shown in Fig. 11 to enable a growth-rate 
comparison. The disturbance is slightly more destabi- 
lized in the full 3-D flow than in the 2-D flow approxi- 
mation. Similar qualitative differences are evident when 
disturbance growth rates in quasi-parallel flow are com- 
pared with those in nonparallel flows. Finally, normal- 
ized primary wave velocity profiles are compared in Fig. 
12. The shapes of the compared profiles agree well. The 
results demonstrate that 2-D simulations capture the 
qualitative features of the true 3-D flow; in addition, be- 
cause a third degree of freedom (w, z) is not present in 
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Fig. 12. Comparison of normalized 2-D and 3-D pri- 
mary wave velocity profiles at A' = 100 in attachment- 
line boundary layer at R = 570 and u> = 0.1249. 
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In the nonparallel theory of Hall et a.l. 8 , the z- 
dependent. form for the flow-accelerated velocity com- 
ponent. w was a key assumption, which led to a system 
of ordinary differential equations rather than partial dif- 
ferential equations. Figure 13 shows the maximum am- 
plitudes of the flow-accelerated velocity component at 
A' = 100 with distance from the attachment line. For 
the present harmonic source, this r-dependent. distur- 
bance form assumed by Hall et a.l. 8 is realized in the 
simulation near the attachment line; however, because 
the harmonic source has a finite length, the disturbance 
behavior near the harmonic-source ends deviates from 
the assumed z dependence. The harmonic-source ends 
cause a perturbation to the flow. Similar difficulties in 
disturbance initialization can be found in the experi- 
ments; however, the core of the test region (i.e. , the at- 
tachment. line) is not. significantly contaminated by t.h|fje 
end effects. 
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Fig. 13. Maximum flow-a.ccelera.t.ed disturbance veloc- 
ity w with distance from attachment, line at A' = 100, 
R= 570, and w = 0.1249. 



In the next, section, the simulation results are pre- 
sented for the neut.ra.l-curve regions described by the 
nonparallel theory of Hall et. a.l. 8 


5.2 The Neutral-Curve Region 


In parameter regimes near the neutral curve, finite 
Reynolds number instabilities are studied near the up- 
per branch, the lower branch, and the critical point, to 
locate the neutral curve. The harmonic-source distur- 


bance generator is used to generate the 2-D modes on 
the attachment, line. 


For the Reynolds number R = 684.2, the 3-D sim- 
ulations yield an upper branch mode for the frequency 
u> = 0.1263; the nonparallel theory of Hall et. a.l. 8 and 
the 2-D DNS results indicate that the upper branch is 


between u> = 0.1230 and u> = 0.1240. Although the 2-D 
and 3-D results yield different, upper branch locations, 
the relative error, or difference, in the locations is only 
about. 2 percent.. This difference may be attributed to 
the assumption that a. 2-D disturbance is generated from 
a. 3-D harmonic source or that, the 3-D base flow does 
not. support, pure 2-D disturbances. 

Near the critical-point, region, digitized data, from 
the results of Hall et. a.l. 8 indicate that the Reynolds 
number R = 580 and frequency u> = 0.1104 is the point, 
furthest, upstream at which an infinitesimal 2-D distur- 
bance becomes unstable. Although this value is not. the 
exact, critical point., this Reynolds number and frequency 
combination lies on the neutral curve in the region of the 
critical point.. The 3-D simulation results suggest, that 
for the frequency u> = 0.1104, the Reynolds number for 
neutral stability is slightly greater than R = 585. This 
leads to less than a. 1-percent, difference between the non- 
pa.ra.llel theory and the simulation results. 

Finally, simulations are performed near the lower 
branch of the neutral curve. The results indicate that 
for the Reynolds number R = 684.2 the lower branch 
of the neutral curve is approximately at the frequency 
l o = 0.082, which agrees with nonparallel theory. 

For practical engineering purposes, the nonparallel 
theory of Hall et. a.l. 8 agrees with the 3-D DNS results 
in the limit, of infinitesimal qua.si-2-D disturbances that 
propagate along the attachment, line. 

5.3 Three-Dimensional Disturbances 

To generate 3-D disturbances, the flow-acceleration 
length of the harmonic-source generator is reduced to 
enable a. more direct, transfer of energy to the w velocity 
component.. Disturbances computed in the parameter 
regime described by a. Reynolds number R = 570 and 
frequency u> = 0.1249 are shown in Fig. 14. By re- 
ducing the length of the original harmonic source from 
-44.2 < A < 44.2 to -20.4 < A < 20.4, the gener- 
ated disturbance is very similar to the previous qua.si- 
2-D disturbance. However, by reducing the harmonic- 
source length to —13.4 < A < 13.4 (one-third of the 
original length), the generated disturbance is signifi- 
cantly stabilized on the attachment, line. The evolu- 
tion no longer represents a. qua.si-2-D disturbance and 
becomes more comparable to a. harmonic point, source. 
Two-dimensional instabilities are apparently dominant, 
on the attachment, line. 
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Fig. 14. Evolution of disturbances in attachment-line 
boundary layer at R = 570 and u> = 0.1249, where dis- 
turbances are generated with harmonic sources of vari- 
ous lengths. 
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Fig. 15. Evolution of disturbance velocity u on at- 
tachment line and top view of 3-D traveling wave in 
attachment-line boundary layer at R = 570 and u> = 
0.1249. 


Next, a harmonic-source generator is used to intro- 
duce a disturbance off the attachment line to determine 
the direction and rate of instability growth or decay. 
The results of a disturbance generated with a harmonic 
source located at 27.8 < Z < 0 are shown in Fig. 15. 
The top view indicates that the harmonic source gener- 
ates a disturbance that evolves along i lie attachment 
line with spreading both away from and toward th§ 
attachment line. These results suggest that the flow- 
accelerated shear away from the attachment line has in- 
sufficient strength to deter the spreading of the distur- 
bance toward the attachment line. Figure 15 also shows 
that the maximum-amplitude u velocity on the attach- 
ment, line initially undergoes a slight decay and then 
continues to grow. The amplitude informattion along 
the attachment line alone indicates an unstable mode 
in the simulations; however, the top view of the flow 
Held indicates that this amplification is caused by the 
wave-spreading phenomenon. The combined amplitude 
and visual results imply that a disturbance generated off 
(but near) the attachment line can supply energy to the 
attachment region by the spreading of the wave pattern. 
In turn, this energy supply may feed an unstable mode 
on the attachment line. 


For the final simulation in this subsection, the 
Reynolds number R = 684.2 and the frequdsjcy u> = 
0.1150 are used because nonparallel theory predicts that 
infinitesimal 2-D disturbances are unstable on the at- 
tachment, line. The disturbance is generated with a, har- 
monic source that, is positioned at, 35.6 < Z < 6.6 (i.e. , 
completely off the attachment, line). The top view of 
the computed disturbance is shown in Fig. 16. The har- 
monic source has generated a, disturbance with a, circu- 
lar pattern. As before, the disturbance evolves primarily 
along the attachment, line, and the wave spreads both 
away from and toward the attachment, line. The ampli- 
tude of the disturbance at, various Z locations is shown 
in Fig. 17. The disturbance has a, peak amplitude at, ini- 
tiated at, Z = 20.4 and undergoes a, strong decay along 
the attachment, line, although the mode is predicted to 
be unstable on the attachment, line. The spread of the 
disturbance toward the attachment, line indicates that, 
the disturbance on the attachment, line is either unsta- 
ble or merely gaining energy at, a, rate comparable to the 
spreading rate. However, because the theory for 2-D dis- 
turbances indicates that, the disturbance is unstable on 
the attachment, line, some combination of energy trans- 
fer due to spreading and linear growth is likely. How- 
ever, the more stable 3-D modes may rob the 2-D mode 
of energy and prevent, flow transition along the attach- 
ment, line. Note that, the u velocity components at, all Z 
locations indicate increased amplitudes with along the 
attachment, line, except, for the Z = 20.4 location which 
indicates decay. Spreading causes the other Z locations 
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to receive energy, but because the Z = 20.4 location was 
the location of maximum initial amplitude and because 
the disturbance propagates along and away from the at- 
tachment line, the location of the maximum velocity is 
no longer at Z = 20.4. This results in an observed de- 
cay at the Z = 20.4 station. Figure 18 shows velocity 
profiles at various Z locations at A' = 100. As energy is 
transferred because of this spreading, the profiles near 
the attachment line undergo a distortion near the wall. 
This distortion leads to multiple maximums and profile 
shapes that deviate from the linear theory. 
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Fig. 18. Comparison of 3-D disturbance velocity profiles 
at A' = 100 near attachment line at R = 684.2 and 
w = 0.1150. 
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Fig. 16. Top view of disturbance evolution in 
attachment-line boundary layer at R = 684.2 and u> = 
0.1150, where disturbance is generated with harmonic 
source near attachment line. 
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Fig. 17. Evolution of disturbance generated off attach- 
ment line in attachment-line boundary layer at R = 
684.2 and w = 0.1150. 


5.4 Suction and Blowing Effects 

By changing the boundary conditions in equation 
(7) from k = 0, steady suction (k < 0) or blowing 
(k > 0) can be used to alter the growth or decay 
of disturbances in the attachment-line boundary-layer 
flow. Near the upper branch of the neutral curve, the 
Reynolds number R = 684.2 and frequency u> = 0.1230 
are used for the simple test case of linear stability with 
suction and blowing. The results of the quasi-2-D dis- 
turbance generated with the elongated harmonic source 
(—44.2 < Z < 44.2) indicate that suction stabilizes and 
blowing destabilizes the disturbance, which agrees with 
the theoretical results by Hall et. al. 8 and the 2-D DNS 
results. 

The results for the 3-D disturbance generated with 
a harmonic source of length 35.6 < Z < 6.6 at the 
Reynolds number R = 684.2 and frequency u> = 0.1150 
indicated growth in the energy on the attachment line 
(Figs. 16 and 17). Because 2-D disturbances that cor- 
respond to this Reynolds number and frequency are lin- 
early unstable on the attachment line, the presence of 
energy should lead to instability growth. Computations 
with suction are used to stabilize the disturbance on and 
near the attachment line. Clearly, Fig. 19 shows that 
suction stabilizes the disturbances located both on and 
off the attachment line. 
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growing disturbances. 

The neutral-curve location predicted by the non- 
parallel theory of Hall -et. al. 8 agreed well with the 3-D 
simulation results in the limit of the infinitesimal 2-D 
disturbances, that propagate along the attachment line. 

For the parameter regions studied here, instabili- 
ties that are generated from harmonic sources located 
off the attachment line spread both toward and away 
from the attachment line. Because of this spreading, 
energy from the initial disturbance is transferred to the 
attachment-line instabilities; however, suction stabilizes 
these instabilities. Furthermore, 3-D instabilities were 
more stable than 2-D or quasi-2-D instabilities. 
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Fig. 19. Effect, of suction on evolution of disturbance 
generated off attachment, line in attachment-line bound- 
ary layer at R = 684.2 and u> = 0.1150. 

6. Concluding Remarks 

In this paper, the results of two- (2-D) and three- 
dimensional (3-D) spatial direct, numerical simulations 
(DNS) of attachment-line instabilities in swept. Hiemenz 
flow are presented. With a. qua.si-para.llel base flow 
approximation, the small- amplitude disturbances were 
shown to grow and decay in agreement, with linear sta- 
bility theory. The true swept. Hiemenz base flow leads to 
a. destabilization of the disturbances, which agrees with 
the nonpara.llel theory of Hall et al. 8 Furthermore, the 
effect, of steady suction and blowing on sma.ll-a.mplit.ude 
disturbances was evaluated with DNS. In agreement, 
with the results of Hall et. al., 8 suction stabilizes and 
blowing destabilizes sma.ll-a.mplit.ude disturbances. 

Subcritica.l instability growth was detected with 
the weakly nonlinear theory of Hall and Malik 9 and 
later refuted with 2-D DNS results by Jimenez et. al. 14 
Theofilis. 15 Although it. is clear that Theofilis 15 did not. 
obtain subcritica.l growth because the initial disturbance 
amplitudes were too small, it. is not. clear why Jimenez 
et. al. 14 failed to obtain subcritica.l growth. The present, 
simulations have found subcritica.l instability growth 
and confirm the former theoretical results of Hall and 
Malik 9 . We surmise that, (similar to nonpara.llel effects) 
nonlinear disturbances broaden the neutral curve to- 
ward higher frequencies in the critical region. Further- 
more, the DNS results demonstrate that, small amounts 
of steady suction stabilize the otherwise sub critically 
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